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Abstract 

We consider a basic model of a dynamical distribution network, modeled as 
a directed graph with storage variables corresponding to every vertex and 
flow inputs corresponding to every edge, subject to unknown but constant 
inflows and outflows. As a preparatory result it is shown how a distributed 
proportional-integral controller structure, associating with every edge of the 
graph a controller state, will regulate the state variables of the vertices, ir- 
respective of the unknown constant inflows and outflows, in the sense that 
the storage variables converge to the same value (load balancing or consen- 
sus). This will be proved by identifying the closed-loop system as a port- 
Hamiltonian system, and modifying the Hamiltonian function into a Lya- 
punov function, dependent on the value of the vector of constant inflows and 
outflows. In the main part of the paper the same problem will be addressed 
for the case that the input flow variables are constrained to take value in an 
interval. We will derive sufficient and necessary conditions for load balanc- 
ing, which only depend on the structure of the network in relation with the 
flow constraints. 
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1. Introduction 

In this paper we study a basic model for the dynamics of a distribution 
network. Identifying the network with a directed graph we associate with 
every vertex of the graph a state variable corresponding to storage, and 
with every edge a control input variable corresponding to flow, which is 
constrained to take value in a given closed interval. Furthermore, some of 
the vertices serve as terminals where an unknown but constant flow may enter 
or leave the network in such a way that the total sum of inflows and outflows 
is equal to zero. The control problem to be studied is to derive necessary 
and sufficient conditions for a distributed control structure (the control input 
corresponding to a given edge only depending on the difference of the state 
variables of the adjacent vertices) which will ensure that the state variables 
associated to all vertices will converge to the same value equal to the average 
of the initial condition, irrespective of the values of the constant unknown 
inflows and outflows. 

The structure of the paper is as follows. Some preliminaries and nota- 
tions will be given in Section 2. In Section 3 we will show how in the absence 
of constraints on the flow input variables a distributed proportional-integral 
(PI) controller structure, associating with every edge of the graph a controller 
state, will solve the problem if and only if the graph is weakly connected. This 
will be shown by identifying the closed-loop system as a port-Hamiltonian 
system, with state variables associated both to the vertices and the edges of 
the graph, in line with the general definition of port-Hamiltonian systems on 

graphs 0, a a a 

; see also jsl, The proof of asymptotic load balancing 
will be given by modifying, depending on the vector of constant infiows and 
outfiows, the total Hamiltonian function into a Lyapunov function. In ex- 
amples the obtained Pl-controller often has a clear physical interpretation, 
emulating the physical action of adding energy storage and damping to the 
edges. 

The main contribution of the paper resides in Sections 4 and 5, where the 
same problem is addressed for the case of constraints on the fiow input vari- 
ables. In Section 4 it will be shown that in the case of zero infiow and outfiow 
the state variables of the vertices converge to the same value if and only if 
the network is strongly connected. This will be shown by constructing a 
Lyapunov function based on the total Hamiltonian and the constraint values. 
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This same construction will be extended in Section 5 to the case of non-zero 
inflows and outflows, leading to the result that in this case asymptotic load 
balancing is reached if and only the graph is not only strongly connected but 
also balanced. Finally, Section 6 contains the conclusions. 

Some preliminary results, in particular concerning Section 3, have been 
already reported before in [7|. 



2. Preliminaries and notations 

First we recall some standard definitions regarding directed graphs, as can 
be found e.g. in (8|. A directed graph Q consists of a finite set V of vertices 
and a finite set £ of edges, together with a mapping from £ to the set of 
ordered pairs of V, where no self-loops are allowed. Thus to any edge e & £ 
there corresponds an ordered pair (f ,w) G V x V (with v ^ w), representing 
the tail vertex v and the head vertex w of this edge. 

A directed graph is completely specified by its incidence matrix B, which 
is an X m matrix, n being the number of vertices and m being the number of 
edges, with (^, j)*'^ element equal to 1 if the j^^ edge is towards vertex i, and 
equal to —1 if the j*^ edge is originating from vertex i, and otherwise. Since 
we will only consider directed graphs in this paper 'graph' will throughout 
mean 'directed graph' in the sequel. A directed graph is strongly connected if 
it is possible to reach any vertex starting from any other vertex by traversing 
edges following their directions. A directed graph is called weakly connected 
if it is possible to reach any vertex from every other vertex using the edges 
not taking into account their direction. A graph is weakly connected if and 
only if keri?"^ = spanl„. Here In denotes the n-dimensional vector with all 
elements equal to 1. A graph that is not weakly connected falls apart into 
a number of weakly connected subgraphs, called the connected components. 
The number of connected components is equal to dim ker B^. For each vertex, 
the number of incoming edges is called the in- degree of the vertex and the 
number of outgoing edges its out-degree. A graph is called balanced if and 
only if the in-degree and out-degree of every vertex are equal. A graph is 
balanced if and only if 1„ G keiB. 

Given a graph, we define its vertex space as the vector space of all func- 
tions from V to some linear space TZ. In the rest of this paper we will take for 
simplicity 7^ = M, in which case the vertex space can be identified with M". 
Similarly, we define its edge space as the vector space of all functions from 
£ to TZ = M., which can be identified with M"^. In this way, the incidence 
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matrix B of the graph can be also regarded as the matrix representation of 
a hnear map from the edge space to the vertex space M". 
Notation: For a, 6 G the notation a ^ b will denote elementwise in- 
equality tti < bi, i = 1, . . . ,171. For ai < bi, i = 1, . . . ,m the multidimensional 
saturation function sat (a; ; a, 6) : M™" — )■ is defined as 

if Xi < tti, 

sat(x ; a, b)i = <^ Xi ii ai < Xi < bi, i = 1, . . . , m. (1) 

if Xi > bi, 

3. A dynamic network model with PI controller 

Let us consider the following dynamical system defined on the graph Q 

X = Bu, X eW\ueW' 

(2) 

y = B^%{x), yew-. 



where if : M" — i- M is any differentiable function, and ^{x) denotes the 
column vector of partial derivatives of H. Here the i^^ element Xi of the 
state vector x is the state variable associated to the i*^ vertex, while Uj is 
a flow input variable associated to the j^^ edge of the graph. System ([2]) 



tne 7 

defines a port-Hamiltonian system (jo], [lo[), satisfying the energy-balance 

±H^u-y. (3) 

Note that geometrically its state space is the vertex space, its input space is 
the edge space, while its output space is the dual of the edge space. 

Example 3.1 (Hydraulic network). Consider a hydraulic network, modeled 
as a directed graph with vertices (nodes) corresponding to reservoirs, and 
edges (branches) corresponding to pipes. Let Xi he the stored water at vertex 
i, and uj the flow through edge j . Then the mass-balance of the network is 
summarized in 

X = Bu, (4) 

where B is the incidence matrix of the graph. Let furthermore H{x) denote 
the stored energy in the reservoirs (e.g., gravitational energy). Then Pi := 
^{x),i = l,...,n, are the pressures at the vertices, and the output vector 
y = B^^{x) is the vector whose j^'^ element is the pressure difference Pi — Pk 
across the j*'* edge linking vertex k to vertex i. 
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As a next step we will extend the dynamical system ([2]) with a vector d 
of inflows and outflows 



X = Bu + Ed, X eW,u eW, de 

(5) 

y = B^^ix), yew-, 



with E an n X k matrix whose columns consist of exactly one entry equal 
to 1 (inflow) or —1 (outflow), while the rest of the elements is zero. Thus E 
specifies the k terminal vertices where flows can enter or leave the network. 

In this paper we will regard as a vector of constant disturbances, and 
we want to investigate control schemes which ensure asymptotic load balanc- 
ing of the state vector x irrespective of the (unknown) disturbance d. The 
simplest control possibility is to apply a proportional output feedback 

OH 

u = -Ry = -RB^ — {x), (6) 

where i? is a diagonal matrix with strictly positive diagonal elements ri, . . . , 
Note that this defines a decentralized control scheme if H is of the form 
H{x) = Hi{xi) + . . . + Hn{xn), in which case the i^^ input is given as times 
the difference of the component of ^{x) corresponding to the head vertex 
of the i^^ edge and the component of ^(x) corresponding to its tail vertex. 
This control scheme leads to the closed-loop system 

f)J-f 

X = ~BRB^ — ix) + Ed. (7) 

ox 

In case of zero in/outflows = this implies the energy-balance 

I.H^-?Ji(.)BRB^?§(.)<,. (8) 

Hence if H is radially unbounded it follows that the system trajectories of 
the closed-loop system ([7]) will converge to the set 

rlT-f 

£ := {X I B^ — {x) = 0}. (9) 

and thus to the load balancing set 

dH 

S = {x\ -T—ix) =al,ae R}. 
ox 



5 



if and only if ker = span{l}, or equivalently |8[, if and only if the graph 
is weakly connected. 

In particular, for the standard Hamiltonian H{x) = this means 

that the state variables Xi{t),i = 1, . . . , n, converge to a common value a as 
t — J- oo. Since ^l^x(t) = it follows that this common value is given as 

For d ^ proportional control u = —Ry will not be sufficient to reach load 
balancing, since the disturbance d can only be attenuated at the expense of 
increasing the gains in the matrix R. Hence we consider proportional-integral 
(PI) control given by the dynamic output feedback 



Xq 

u 



-Ry - ^(^c) 



(10) 



where Hc{xc) denotes the storage function (energy) of the controller. Note 
that this PI controller is of the same decentralized nature as the static output 
feedback u = —Ry. 

The j^^ element of the controller state Xc can be regarded as an additional 
state variable corresponding to the j*^ edge. Thus Xc G M™, the edge space 
of the network. The closed-loop system resulting from the PI control ( ITOl) is 
given as 

" -5 



X 




Xc 





-BRB^ 
B^ 







dx 

dHc^ 

. dxc 



X) 



+ 



d. 



(11) 



This is again a port-Hamiltonian systeno, with total Hamiltonian Htot{x, Xc) 
:= H{x) + Hc{xc), and satisfying the energy-balance 



d 
It 



tot 



;i2) 



Consider now a constant disturbance d for which there exists a matching 
controller state Xc, i.e., 

Ed = B^(xc). (13) 



dx 



3 See for a general definition of port-Haniiltonian systems on graphs. 
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This allows us to modify the total Hamiltonian Htot{x,Xc) intc0 

Vd{x,Xc) := H{x) + Hc{xc) 7: — -{xc){xc - Xc) - Hc{xc), (14) 

OXc 

which will serve as a candidate Lyapunov function; leading to the following 
theorem. 

Theorem 1. Consider the system ^ on the graph Q in closed loop with the 
Pl-controller ^^). Let the constant disturbance d he such that there exists a 
Xc satisfying the matching equation l[T3\) . Assume that radially 
unbounded. Then the trajectories of the closed-loop system ( fJTl) will converge 
to an element of the load balancing set 

dH dH - 

Stot = {{x, Xc) I — (x) = al, a e R, B—^{xc) = Ed }. (15) 

ox OXr 



if and only if Q is weakly connected. 

Proof. Suppose Q is weakly connected. By f|T2|) for d = d we obtain, making 
use of ([I3D, 

rf,, d'^H^ ,^^^rdH, , H , 
-V,= ~[.)BRB--(.) + —(.)Ed- 

^(R)B'-f^(x) (16) 

OXc OX 

= -^ix)BRB^^{x)<0. 
ox ox 

Hence by LaSalle's invariance principle the system trajectories converge to 
the largest invariant set contained in 

OH 

{{x,Xc)\B''—{x) = 0}. 

Substitution of B'^^(x) = in the closed-loop system equations (|TT1) yields 
Xc constant and —B^'^^{xc) + Ed = 0. Since the graph is weakly connected 
B^^{x) = implies ^{x) = at. If the graph is not weakly connected then 
the above analysis will hold on every connected component, but the common 
value a will be different for different components. □ 



This function was introduced for passive systems with constant inputs in 
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Corollary 2. If kei B = 0, which is equivalent (f^]) to the graph having no 
cycles, then for every d there exists a unique Xc satisfying / f73]) . and conver- 
gence is towards the set £^tot = {{x, Xc) \ ^{x) = at, a E Xc = Xc}- 

Corollary 3. In case of the standard quadratic Hamiltonians H{ 

Hc{xc) = |||a^c|P there exists for every d a controller state Xc such that /[T^) 

holds if and only if 

imEcimE. (17) 

Furthermore, in this case Vj equals the radially unbounded function + 
|||xc — XcW^, while convergence will be towards the load balancing set Stot = 
{{x,Xc) \ X = al,a eM, Bxc = Ed}. 

A necessary (and in case the graph is weakly connected necessary and 
sufficient) condition for the inclusion imii^ C im B is that I'^E = 0. In its 
turn I'^E = is equivalent to the fact that for every d the total inflow into 
the network equals to the total outflow). The condition I'^E = also implies 

t^x = -t^BRB^^(x) + l^Ed = 0, (18) 
ox 

implying (as in the case d = 0) that I'^x is a conserved quantity for the closed- 
loop system ffTTj) . In particular it follows that the limit value \imt^oox{t) e 
span{l} is determined by the initial condition x(0). 

Example 3.2 (Hydraulic network continued). The proportional part u = 
—Ry of the controller corresponds to adding damping to the dynamics (pro- 
portional to the pressure differences along the edges). The integral part of 
the controller has the interpretation of adding compressibility to the hydraulic 
network dynamics. Using this emulated compressibility, the Pl-controller is 
able to regulate the hydraulic network to a load balancing situation where 
all pressures Pi are equal, irrespective of the constant inflow and outflow 
d satisfying the matching condition [W\) . Note that for the Hamiltonian 
H{x) = |||a;|p the pressures Pi are equal to each other if and only if the 
water levels Xi are equal. 

4. Constrained flows: the case without in/out flows 

In many cases of interest, the elements of the vector of flow inputs u G M"* 
corresponding to the edges of the graph will be constrained, that is 

u G W := {m G I ^ u < M+} (19) 
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for certain vectors u~ and satisfying u~ ^ ^ uf ,i = 1, . . . ,m (through- 
out ^ denotes element- wise inequality). This leads to the following con- 
strained versioijf] of the PI controller ffTOl) given in the previous section 

= y, 

(20) 

u = sat [- Ry - ^{xc) ■,u ,u+) 

Throughout this paper we make the following assumption on the flow con- 
straints. 

Assumption 4. 

ur ^ 0, uf ^0, <ut,i = l,...,m (21) 

It is important to note that we may change the orientation of some 
of the edges of the graph at will; replacing the corresponding columns bi 
of the incidence matrix B by —bi. Noting the identity sa.t{—x;u~,u^) = 
— sat(2; ; —uf, —u~) this implies that we may assume without loss of gener- 
ality that the orientation of the graph is chosen such that 

u-^0<uf, t = l,...,m (22) 

This will be assumed throughout the rest of the paper. In general, 
we will say that any orientation of the graph is compatible with the flow 
constraints if (122|) holds. If the j-th edge is such that uj = then we will 
call this edge an uni- directional edge, while if m~ < then the edge is called 
a bi-directional edge. 

In this section we will first analyze the closed-loop system for the con- 
strained Pl-controller under the simplifying assumption of zero inflow and 
outflow {d = 0). In the next section we will deal with the general case. 
Furthermore, for simplicity of exposition we consider throughout the rest of 
this paper the standard Hamiltonian Hc{xc) = llja^dP for the constrained 
PI controller and the identity gain matrix R = I, while we also throughout 
assume that the Hessian matrix of Hamiltonian H{x) is positive definite for 
any x. Thus we consider the closed-loop system 

X = Bsat{-B^^{x)-x,;u-,u+), 

X, = E^f (a:). ^^^^ 



^See also [12] for a related problem setting where a constrained version of the propor- 
tional controller ([6|) is considered. 
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In order to state the main theorem of this section we need one more definition 
concerning strong connectedness with respect to flow constraints. 



Definition 5. Consider the directed graph Q together with the constraint 
values u~,u~^ satisfying ^ME)- Then we will call the graph strongly connected 
with respect to the flow constraints u~ ^ u ^ if the following holds: for 
every two vertices vi,V2 there exists an orientation of the graph compatible 
with the flow constraint^ and a directed path (directed with respect to this 
orientation) from vi to V2- 



Theorem 6. Consider the closed-loop system fl23|) on a graph Q with flow 
constraints ^ u ^ satisfying ^W) . Then its solutions converge to the 
load balancing set 

OH 

£tot = {ix,Xc) I -t:—{x) = aln, B sat(-x c;u~,u^) = 0} (24) 
ox 

if and only if the graph is strongly connected with respect to the flow con- 
straints. 

Proof. Sufficiency: Consider the Lyapunov function given by 

OH 

V{x, X,) = tlS{- B^-^{x) - Xe ; M", u+) + H{x), (25) 

with 

S{x;u',u'^)i := sat{y;Ui,uf)dy. (26) 



It can be easily verifled that V is positive deflnitive, radially unbounded and 
C^. Its time-derivative is given as 

= — sat ( — -B (a;) — aJc; M~, M j-B -Bsat( — -B — — (x) — Xc ; m", m j 
^ 0. 

(27) 



^Note that for different pairs of vertices we may need different orientations compatible 
witli tlie flow constraints. Thus the definition of strong connectedness with respect to the 
flow constraints is weaker than the existence of an orientation of the graph compatible 
with the flow constraints in which the graph is strongly connected. 
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By LaSalle's invariance principle, all trajectories will converge to the largest 
invariant set, denoted as X, contained in /C = {{x,Xc) \ Bsat ( — B'^^{x) — 
"'") = 0}. Whenever x G /C it follows that x = and thus x{t) = v 
for some constant vector v. Hence, since Xn = B^^(x), it follows that 

x,{t) = B^^iu)t + xM- 



dx 

dx 



Suppose now that B^^{iy) ^ 0. Then for t large enougn 



= sat ( - B^^iu) - B^^{u)t - xM,u-,u^) 

ox ox ox 

- 1 (^^f ('')).^.. 



(28) 



where 



ut 



if(5-f(.)):<0. (29) 



, m. 



Hence , in view of ^ < u^, we have (i?"^^(z/)) . ^ 0, for i = I, 
However since the graph is strongly connected with respect to the flow con- 
straints, if {B^^{iy)). > 0, then there exists j such that {B^^{p)). < 0. 

This yields a contradiction. We conclude that B^^(iy) = 0, which implies 
^(z/) = atn, and thus all trajectories converge to £^tot- 

Necessity: Assume without loss of generality that the graph is weakly 
connected. (Otherwise the same analysis can be performed on every con- 
nected component.) If the graph is not strongly connected with respect to 
the flow constraints then there is a pair of vertices Vi, Vj for which there exists 
a compatible orientation and a directed path from Vi to vj, but not a com- 
patible orientation and directed path from vj to Vi. In other words, there can 
be positive flow from Vi to Vj, but not vice versa. Then for suitable initial 
condition, |^(x(t) < ^{x{t) for all t ^ 0, and thus there is no convergence 

to Stot- □ 



Remark 7. Note that for — )■ — oo, 6j — )■ oo, the Lyapunov function ( 125|1 
tends to the function H(x) + ^\\B^^(x) +Xc|p, which is different from the 
Lyapunov function H{x) + |||xc|p used in the previous section. 

In the special case that the flow constraints are such that all the flows Ui 
can follow both directions, we obtain the following corollary. 

Corollary 8. For a network with constraint intervals [u'[, uf] with < < 
uf , i = 1, . . . ,m, the trajectories of the closed-loop system f l23|) will converge 
to the set if and only if the network is weakly connected. 
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Proof. In this case (since all the edges are bi-directional) weak connectedness 
is equivalent to strong connectedness with respect to the flow constraints. If 
the graph is not weakly connected then the components of ^ will only 
converge to a common value on every connected component. □ 



5. Nonzero inflows and outflows 

In this section we deal with the general case of nonzero (but constant) 
inflows and outflows d. Thus we consider the closed-loop system 

dH 

X = BsaX { — B'^——(x) — Xc ; u~, u"*") + Ed, 

(30) 



Xr. = B^^ix) 



dx 



with im.Ec im.B. 

In order for the system to reach consensus, we need to impose conditions 
on the magnitude of the in/outflows d. 

Definition 9. Given the constraint values u~ < the permission set 
V{u^, u^) are defined as 

where the intervals Vi{u~,u~^) is defined by: 

( {u7,-ui) ifOe {ui,ui) and ^ 
Pi(M~,u+) = <j {-ut,ut) ifOe{Ui,uf)and\u:r\>\uf\ (31) 
( (0, u+iJ if (u", M+) = (0, uf), 

where = minjw^ | i such that = 0}. 

Theorem 10. Consider a graph Q with dynamics Suppose that every 

edge allows bi-directional flow, i.e., u~ < < uf, i = 1, . . . ,m. Then for any 
in/outflow d for which there exists x^ G V{u~,u~^) such that Ed = Bxc, the 
trajectories of (E^i ''^M converge to 



dH 

S^^i = {{x,Xc) I = al,ae R, B sat{-Xc]u~ ,u+) + Ed = 0}. (32) 

if and only if the graph Q is weakly connected. 
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Proof. By the matching condition Ed = Bxc and the identity 



sat(x — ?7;-u , -u^) + = sat(x ; m + r], + r^), Vr] G M", (33) 
the system (15U]1 can be written as 



dH 

X = B Sat(— 5"^-^— (x) — Xc'jU^ + Xc, + Xc) 

Xc = B -7—[x), 
ox 



(34) 



where Xc = Xc — Xc- Since by construction {u + Xc)i < 0, i = 1, . . . , m, the 
proof now follows from Corollary [SI □ 

The following theorem covers the case that every edge is uni-directional. 

Theorem 11. Consider a network Q with dynamics / TfOj) with flow con- 
straints such that = 0,i = 1, . . . ,m (uni-directional flow). Then for any 
G M™ and any in/outflow d for which there exists Xc G V{Om,u~^) such 
that Ed = Bxc, the trajectories of [3U\} converge to 



dH 

^tot = {(x, Xc) I -T^ix) = at, a E M, B sat(— Xc ; 0^, m^) + Ed = }, (35) 
ox 

if and only if the graph in the (only) orientation compatible with the flow 
constraints is strongly connected and balanced. 

In order to prove Theorem [TT] we need the following lemma. Recall that a 
directed graph is balanced if every vertex has in-degree (number of incoming 
edges) equal to out-degree (number of outgoing edges). Furthermore, we will 
say that two cycles of a graph are non- overlapping if they do not have any 
edges in common. 

Lemma 12. A strongly connected graph is balanced if and only if it can be 
covered by non- overlapping cycles. 

Proof. Sufficiency: If a graph can be covered by non-overlapping cycles, then 
every vertex necessarily has the same in-degree and out-degree; so this graph 
is balanced. 

Necessity: Since the graph is strongly connected, every two vertices can 
be connected by a directed path, and the graph can be covered by cycles. 
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Now suppose that the graph cannot be covered by non-overlapping cycles. 
We will show that this implies that the graph is not balanced. 

Let k be the smallest number of cycles needed to cover the graph, and let 
T = (Ci, C2, . . . , Ck) be a covering set of cycles. According to our assump- 
tion, at least one edge of the graph is shared by two or more cycles in T. 
We claim that the set of shared edges can not contain any cycles. Indeed, 
suppose there is one cycle, denoted as V (depicted in FigdJ^a)), whose edges 
are all shared by elements of T. If P G T, then obviously T is not a minimal 
covering set, since by deleting the cycle V from T we have a covering set of 
k — 1 elements. 

Thus T> ^ T. It can be seen that the minimal number c of cycles in T 
which cover T> twice is at least 4. Denote such a minimal set of c cycles in 
T which cover P by Tp. We will now show that by combining these c cycles 
with the cycle V there exist 3 cycles in the original graph Q which cover the 
subgraph given by Tv', thus reaching a contradiction with the minimality of 
T. The construction of these 3 cycles is indicated in Figure [1] Consider 
for simplicity the case that 4 cycles in T, denoted by Ci, C2, C3, C4 cover V 
twice. Combining (depending on the orientation of the cycles) part of Ci 
with part of C3, and part of C2 with part of C4 (see Figured]), we can define 
2 cycles which together with the cycle D yields a set of 3 cycles which cover 
the subgraph spanned by Ci, C2, C3, C4. 

In conclusion, there must exist at least one shared edge, say (fj,fj), such 
that all edges with tail-vertex Vj are used only once in T. But this implies 
that Vj has larger out-degree than in-degree, i.e., the graph is unbalanced. □ 



Proof of Theorem [TTl Sufficiency: By using (133|1 we rewrite the system as 

(36) 



i; = i? sat ( — 5"^^— (x) — Xc ; Xc, u'^ + Xc) , 



where Consider now the Lyapunov function 

rlT-f 

V{x, Xc) = t^S{- B^—{x) - Xc ; Xc, u+ + Xc) + H{x). (37) 

Similar to the proof of Theorem El if a solution {x{t),Xc(t)) is in the largest 
invariant set X contained in {{x,Xc) 1^ = 0}, then x is a constant vector. 
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(a) The cycle V split 
into two parts Pi, I?2- 



(b) The subgraph To 



Figure 1: (a). The cycle T). We divide the edges of V into two disjoint sets: "Di contains 

the left part of T) and 1^2 contains the rest. (b). The subgraph given by Tv- Without T>, 
we need at least 4 cycles to cover V twice; these cycles are given as Ci = Tvi U 2?i,C2 = 
Tdt. U X'ljCs = T-D3 U V2 and C4 = Tda U I'2- It follows that To is also covered by the 3 
cycles: ©(clockwise), 7x)2 U Tea (counterclockwise) and 7x>i U7i)4 (counterclockwise). 

denoted as v. Furthermore, I is given as 



X = {(i/, X,) I X, = B^—{v)t + Xe(0), 
Ssat ( - B^—{u) - B^—{v)t - Xc(0) ; Xc, 1*+ + Xc) = 0, > 0}. 




(39) 



j=i 



where 




(40) 



Since the graph is balanced we have B\rn — 



0, and thus 




(41) 
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E 



By the definition of the permission set V{Om, W^), < Xd < + Xcj for any 
= 1,2, . . . ,m, so 

{B'^i-))fi<^- (42) 

i=l 

This yields a contradiction. Hence B^^(u) = and therefore 

dH 

X= {{u,Xc)\-^{i^) = cl„,5sat(-Xc;Xc,M+ + = 0} 

Necessity: First, if the graph Q is not strongly connected, by the same argu- 
ment in Theorem |6l it can be easily seen that ^ will not reach consensus. 

Now we will show that if the strongly connected network is unbalanced, 
then there exists a constraint interval [0^, u^] and an in/outflow d for which 
there exists Xc G 7^(0^,^^) such that Ed = Bxc while ^{x) is not converg- 
ing to consensus. 

For simplicity of exposition we shall take the constraint interval as [Om, Im]- 
As in the proof of Lemma [T2] we let k be the minimal number of cycles to 

cover Q, and we let T = (Ci, . . . , Ck) be a minimal covering set for Q. With 

some abuse of notation 

BQ = 0,i = l,...,k (43) 

where Ci is the m-dimensional vector whose j-th component is equal to the 
number of times the j-th edge appears in the cycle Cj. 

In the following, we will prove that there exist B^^{u) ^ 0, Xc G 
V{Qm, Im), Xc{0) and A e M, such that 

OH OH 

sat ( - B^—{v) - B^-Q^{^)t - ; a^c, tm + x,) = AT, Vt > 0, (44) 

where u is the equilibrium value of x as above and T is the m-dimensional 
vector whose i-th component is the number of cycles in T which contains 
i-th edge. This implies that the system has an equilibrium which satisfies 
If (z/) ^ span{l J. 

Consider as above a minimal covering set T = (Ci, . . . , Ck) for Q. Let 
Tmax := {Ti \ i = 1,2, ...,m}, and denote Si = {i-th edge | Tj = Tmax}- 
Every cycle in T has at least one non-overlapped edge (see the proof of 
Lemma [12]), and we denote by S2 the set of all the non-overlapped edges 
in the cycles in T which contain at least one edge which is overlapped /max 
times. 
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dH , 


>; 




< 


dH 


dx 


h 


dx 


dH , 


y, 




> 


dH 


dx 


h 


dx 


dH . 


>; 






dH 


dx ' 


h 




dx 



In the last step, we will make the flows through the edges in 8i reach the 
upper bounds of the constraint intervals, and the flows through the edges in 
S2 reach their lower bounds. By taking 



i, {v„Vj)e£2 (45) 
i, else 

for suitable Xc and Xc(0), it follows that (jS]) holds. Indeed, in the set S1US2, 
the equation ( H^ takes the form 

1 + Xcq = XTq, g — th edge belongs to Si 
Xcp = XTp,p — th edge belongs to £2 

Now take A be such that -r— < A < 1. Then (USD contains |£i| + |£^2l equations 
and the same number of variables, and has a unique solution such that 

< Xcp < 1 , , 

' 47 

0<Xe, <1 ^ ' 

Furthermore, pick Xc{0) in the third equation of ( l45i) as 

x^{0)r = -AT,, r - th edge belongs to S\{Si U S2). (48) 
Obviously, there exists < Xcr < 1 such that 

Xcr < —Xc{0)r < 1 + Xcr (49) 

In conclusion, there exists an equilibrium (z/, Xc) that does not satisfy 
B^^{u) = 0, and thus ^{j^) can not reach consensus. □ 

The above constructive proof is illustrated by the following example. 

Example 5.1. Consider a directed graph in Figure^ with dynamics given 
by system fl30|) where H{x) = and [u~,u^] = [07,17] 

x = B sat(-B'^x - Xc ; O7, 17) + Ed, 

(50) 

Xc = B X. 

The purpose of this example is to show there exist in/outflows d satisfy- 
ing the matching condition for which x does not converge to consensus. 
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Figure 2: Network of Example 15.11 




Timet 

Figure 3: The time-evolutions xi{t),X2{t),X3{t),X4{t),X5{t) of the system (|50p . 

By taking Ed = Bxc where Xc = |l7,x(0) = (3,75,1,4)-^ and Xc{0) = 
(1, —1, —1, 1, 1, 1, 1)^, the state x in system (150|) will converge to v with 

= > ^^5 > > v\ cind z/4 < ui as can be seen from the numerical 
simulation in Figure^ 

The same result can be deduced from the following analysis. In Figure^ 
the smallest number of cycles to cover the whole graph is 3; one option being 
(66,67,62,63), (61,62,63), (63,64,65). So BT = where 

T = (1,2, 3, 1,1, 1,1)^ (51) 

In this case Si = {63}, S2 = {^i, 64, 65, Cq, 67}. By setting 1^2 = ^3 > ^5 > ^4: > 
z/i, the flow in 63 reaches the upper bound, while the flows in 61,64,65,65,67 
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reach the lower bounds, i.e. 



sat{-B^u - B^ut - Xc{0) ; x^, x^ + l) = ^T, Vt > 0. (52) 
Thus there exists an equilibrium v satisfying B^u ^ 0. 



6. Conclusions 

We have discussed a basic model of dynamical distribution networks 
where the flows through the edges are generated by distributed PI controllers. 
The resulting system can be naturally modeled as a port-Hamiltonian sys- 
tem, enabling the easy derivation of sufficient and necessary conditions for the 
convergence of the state variables to load balancing (consensus). The main 
part of this paper focusses on the case where flow constraints are present. A 
key ingredient in this analysis is the construction of a Lyapunov function. 
We distinguish between the case that the flow constraints corresponding to 
all the edges allow for bi-directional flow and the case that all the edges only 
allow for uni-directional flow. For both cases we have derived necessary and 
sufficient conditions for asymptotic load balancing based on the structure of 
the graph. 

An obvious open problem is the extension of our results to the general 
case where some of the edges allow bi-directional flow and others only uni- 
directional flow. This is currently under investigation. Many other questions 
can be addressed in this framework. For example, what is happening if the 
in/outflows are not assumed to be constant, but are e.g. periodic functions 



of time; see already [13||. Furthermore, the use of constrained Pl-controUers 



may suggest a fruitful connection to anti-windup control ideas. 
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